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It has been suggested that the nearest-neighbour antiferromagnetic Heisenberg model on the 
Kagome lattice may be a good starting point for understanding the spin-liquid behaviour discovered 
in Herbertsmithite. In this work, we investigate possible quantum spin liquid phases in the presence 
of spin-rotation symmetry-breaking perturbations such as Dzyaloshinskii-Moriya (DM) and Ising 
interactions, as well as second-neighbour antiferromagnetic Heisenberg interactions. Experiments 
suggest that such perturbations are likely to be present in Herbertsmithite. We use the projective 
symmetry group (PSG) analysis within the framework of the slave-fermion construction of quan- 
tum spin liquid phases and systematically classify possible spin liquid phases in the presence of 
perturbations mentioned above. The dynamical spin-structure factor for relevant spin liquid phases 
is computed and the effect of those perturbations are studied. Our calculations reveal dispersive 
features in the spin structure factor embedded in a generally diffuse background due to the ex- 
istence of fractionalized spin- 1/2 excitations called spinons. For two of the previously proposed 
Z2 states, the dispersive features are almost absent, and diffuse scattering dominates over a large 
energy window throughout the Brillouin zone. This resembles the structure factor observed in re- 
cent inelastic neutron scattering experiments on singlet crystals of Herbertsmithite. Furthermore, 
one of the Z 2 states with the spin structure factor with mostly diffuse scattering is gapped, and it 
may be adiabatically connected to the gapped spin liquid state observed in recent Density Matrix 
Renormalization Group calculations for the nearest neighbour antiferromagnetic Heisenberg model. 
The perturbations mentioned above are found to enhance the diffuse nature of the spin structure 
factor and reduce the momentum dependencies of the spin gap. We also calculate the electron spin 
resonance (ESR) absorption spectra that further characterize the role of spin-rotation symmetry 
breaking perturbations in these states. We suggest that the measurement of ESR spectra can shed 
more light into the nature of the ground state in Herbertsmithite. 



I. INTRODUCTION 

Two dimensional spin- 1/2 Mott Insulators in corner 
sharing lattice geometries, like the Kagome lattice, have 
been studied extensively due to their natural tendency to 
suppress magnetic ordering and potential in realizing a 
quantum spin liquid (QSL) ground state. 1 27 While the 
theoretical understanding of QSLs suggests several inter- 
esting features of these highly entangled quantum states 
of matter, such as effective quantum number fraction- 
alization and/or topological order, 1,2,28,29 its conclusive 
experimental realization has remained elusive so far. 

Of the several candidate materials investigated, the 
compound, Herbertsmithite (ZnCu3(OH)eCl2) seems 
particularly promising. 8,9,11 15 > 26 > 30 33 Here, the Cu 2+ 
ions sit on an isotropic Kagome lattice and carry spin- 1/2 
moments. While the Curie- Weiss temperature (6cw~ 
—300 K) suggests a strong and predominantly antiferro- 
magnetic exchange interaction, the spins do not show 
any sign of long-range order or freezing down to 50 
mK. 14,15 Inelastic neutron scattering on powder sam- 
ples shows magnetic scattering for a broad range of 
wavenumbers; 14,30 furthermore, recent inelastic neutron 
scattering on single crystals show that even at tempera- 
tures as low as 1.6 K, a diffuse magnetic scattering con- 
tinuum exists over a large energy window (1.5 — 11 meV) 



in significant portions of the Brillouin zone and no spin 
gap down to at least 1.5 meV. 9 This has raised hopes 
that a low-temperature QSL phase may be realized in 
this material. In a QSL, such a continuum of scattering 
is naturally expected due to the presence of fractional- 
ized spin-1/2 excitations, called spinons, as opposed to 
the conventional spin-1 magnon excitations in a magnet- 
ically ordered systems. 

While Herbertsmithite suffers from inter-site disorder, 
it seems to mostly come in the form of 14% occupation 
of inter-layer non-magnetic Zn 2+ ions by excess magnetic 
Cu 2+ ions, with less than 5% Zn 2+ impurities in the 
Kagome planes. 8 Hence, it has been suggested that this 
material is, to a very good approximation, a realization 
of spin-1/2 nearest neighbour (NN) Heisenberg antifer- 
romagnet on the Kagome lattice . 14 > 32 > 34 > 35 We term this 
the nearest-neighbour Heisenberg Kagome antiferromag- 
net (HKAF). 

The above suggestion has brought in focus the 
theoretical studies on the NN-HKAF. These have a 
long and varied history, and there exist several pro- 
posed candidate ground states. Such states include 
valence-bond solids, 17,18,36 38 gapless 39 41 and gapped 
QSLs. 10,19,25,42,43 Of these, a gapped QSL was found 
recently by Density Matrix Renormalization Group 
(DMRG) calculations. 5,7,27 However, variational Monte 
Carlo (VMC) studies have found a gapless spin- 
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disordered ground state with energy very similar to that 
of the DMRG. 4 These and other studies also seem to sug- 
gest that the strictly NN-HKAF may be close to a quan- 
tum phase transition. Indeed, recent studies find that 
the nature of the ground state can be changed by tun- 
ing the small next-nearest-neighbour (NNN) exchange 
interaction. 6,44 ' 45 

For Herbertsmithite, this means that, due to the po- 
tential proximity to a quantum phase transition, small 
perturbing interactions may play an important role, par- 
ticularly at low energies. 46 Hence, it may be important 
to account for these smaller perturbing energy scales in 
order to understand the experiments on Herbertsmithite. 
Electron spin resonance (ESR) and magnetic susceptibil- 
ity measurements suggest that perturbations break spin- 
rotation symmetry, with both a Dzyaloshinsky-Moriya 
(DM) interaction and an Ising-like easy- axis exchange, 
each of strengths around 10 percent of the NN antiferro- 
magnetic exchange interaction. 8 ' 13,47 ' 48 

In this paper, we study the effect of these further per- 
turbations on various QSL states that potentially offer 
an explanation of the unusual phenomenology of the non- 
magnetic ground state in Herbertsmithite. We contrast 
these states with respect to their signature in spin-spin 
correlations, which is measured in inelastic neutron scat- 
tering measurements, by calculating the dynamical spin- 
structure factor for a host of candidate Z2 and associated 
U(l) spin liquid states that are allowed by the projec- 
tive symmetry group analysis on the Kagome lattice. 10 
In addition, we calculate the ESR absorption spectra 
that provide useful information in systems without spin- 
rotation invariance, and hence can help us to identify the 
nature of the QSL. 49 ' 50 ESR absorption resulting from 
DM interactions has been probed in, for instance, the 
quasi- two-dimensional system CS2CUCI4. 51 In this com- 
pound, the ESR absorption can be interpreted in terms 
of deconfined spin- 1/2 spinons, in contrast with the spin 
chain Cu Benzoate, where bound states dominate the 
ESR spectra. 52 ' 53 

Due to the involved nature of the results presented 
here, before discussing further details, we briefly summa- 
rize first the current situation of the spin liquid physics 
on Kagome antiferromagnets in view of Herbertsmithite, 
and second the main results obtained in this work. 



A. Summary of results: Spin liquid physics in 
Herbertsmithite 

The current theoretical framework of constructing 
QSL states is largely based on slave-particle mean- 
field theories. 28 ' 29 ' 54 58 In particular, we shall focus on 
the slave-fermion mean-field theory where each spin- 
1/2 operator is written in terms of a fermion bilin- 
ear: Si = \f} ot o' ot pfip (also see Eq. 5) with the con- 
straint J2afiafia = 1- These fermions, called spinons, 
carry spin- 1/2 and form the low energy quasi-particles 
whose excitation spectrum, as we shall see below, mainly 



characterizes the low energy response of a QSL state. 
The slave-particle representation introduces gauge redun- 
dancy (see Eq. (14)), and hence the spinons are mini- 
mally coupled to an emergent gauge field. Depending on 
the structure of the gauge group, we can broadly classify 
quantum spin liquids as SU(2), U(l) or Z2 spin liquids. 

In the mean-field theory, the spin-spin interactions 
give rise to quart ic spinon interactions, which are decou- 
pled in the spinon hopping (particle-hole) and pairing 
(particle-particle) channels (see Section IIIB). These to- 
gether constitute a QSL ansatz, and different mean field 
QSL states are characterized by the structure of these 
ansatze. However, the above gauge redundancy suggests 
that the mean- field states are only projective representa- 
tions of physical states. Hence, the same physical state 
may be represented by a set of gauge equivalent mean- 
field ansatze. 28 ' 29 Different mean- field spin liquid ansatze 
form different projective representations of the lattice 
(and time-reversal) symmetries of the Hamiltonian, and 
the spinons also transform accordingly. Thus, a system- 
atic understanding of distinct QSLs requires a careful 
analysis of the projective symmetry group (PSG). 28 ' 29 In 
the presence of both spinon hopping and pairing terms, 
the mean field calculations for the spinon ground state 
can be cast in the form of a BCS type Hamiltonian for 
the spinons (see Section IIIB, Eq. (11)). If the Hami- 
tonian is spin-rotation symmetric, the hopping and the 
pairing channels can be completely separated into a spin 
independent (singlet) part and a spin dependent (triplet) 
part (Section IIIB). In this work we refer to them as sin- 
glet and triplet channels (or terms) respectively. In usual 
studies of models with frustrated antiferromagnetic spin- 
rotation invariant exchanges, the singlet ansatze give rise 
to stable QSL mean field solutions and hence the triplet 
channels are absent. Thus, we obtain symmetric QSLs 
that do not spontaneously break spin-rotation symme- 
try. However, it is known that in the presence of com- 
peting ferromagnetic and antiferromagnetic interactions, 
both singlet and triplet channels can exist (and may give 
rise to exotic three dimensional gapped U(l) spin liquid 
states with topologically protected gapless surface spinon 
modes 59 ' 60 ). 

A PSG analysis 10 of the singlet ansatze on Kagome lat- 
tice shows that there are eight Z2 QSL states with gapped 
or gapless fermionic spinons that can be stabilized by NN 
and NNN hopping and pairing terms for the spinons. All 
these states are found in the neighbourhood of four par- 
ent U(l) states, one of which, with Dirac-like spinon ex- 
citations, was seemingly realized as the ground state of 
the NN-HKAF model in recent VMC studies. 4 Hence, we 
start by looking at these parent QSLs. Since the above 
spin liquids are time reversal symmetric, there exists a 
gauge where the singlet U(l) QSL mean-field ansatze 
can be written employing NN real, singlet hopping terms. 
These spin liquids are conveniently labelled using the fol- 
lowing notation: J7(l)[a, 6], 10 where a (b) denotes the 
magnetic flux of the emergent gauge field through the 
triangular (hexagonal) plaquette of the Kagome lattice. 
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Time reversal symmetry constrains these fluxes to be ei- 
ther or 7T. Accordingly, we have [7(1) [0,0], f7(l)[0,7r], 
f7(l)[7r,7r], and E/(1)[tt,0] QSLs. 10 ' 39 From these four 
U(l) QSLs, the eight Z2 QSL are obtained by tuning in 
appropriate spinon-pairing terms (up to second nearest 
neighbour) allowed by PSG. These pairing terms neces- 
sarily break the gauge structure from U(l) down to Z2, 
hence the name. The nomenclature for the Z2 spin liq- 
uids is derived straightforwardly from their parent U(l) 
states. 10 The added Latin/Greek letter distinguishes be- 
tween two or more Z2 QSL obtained from the same U(l) 
parent. 10 Both the U(l) and Z2 QSLs have characteristic 
spinon band structure. We note that, while this spinon 
dispersion is not a gauge invariant quantity, it is still use- 
ful to understand its general features in order to capture 
various properties of these states like the nature of the 
experimentally observable spin-structure factor. 61 Hence, 
we summarize the general characteristic features of the 
spinon spectra (discussed in detail in Section IV A) in the 
second column of Table I. 

More recent work, however, indicates that the NN- 
HKAF model likely sits close to a quantum phase tran- 
sition between a Z2 spin liquid and a valence-bond- 
solid state upon tuning the second-neighbour interac- 
tion, or a magnetically ordered state upon tuning the 
Dzyaloshinsky-Moriya interaction. 5 7,11,12,16,46 It is this 
proximity to the quantum phase transition that may 
make the system very sensitive to small perturbations 
which then can have a sizable effect. 46 The perturba- 
tions to the HKAF seen in Herbertsmithite are then likely 
to be relevant to a proper description of the material, 
and we endeavour to gain an understanding of how the 
aforementioned symmetric spin liquid states are affected 
under these perturbations, particularly the breaking of 
spin-rotation symmetry induced by DM interaction and 
Ising anisotropics. 

Since the spin-rotation symmetry is broken, the singlet 
and the triplet spinon decoupling channels (both hopping 
and pairing) can no longer be separated, and a more gen- 
eral QSL ansatze incorporating both singlet and triplet 
decomposition, and allowing their intermixing, needs to 
be introduced. 62 While it is possible to use such a general 
four- component spinor representation, as in the case of a 
triplet BCS Hamiltonian, here, we find that since the DM 
vectors and the Ising anisotropics both are perpendicu- 
lar to the Kagome plane, we can use a two-component 
spinor representation. The effect of such spin-rotation 
symmetry breaking terms on the low-energy excitations 
of these spin liquids is also summarized in the Table I 
(third column). Generally, the triplet terms decrease the 
density of states at the spinon Fermi-level in most of the 
17(1) and Z 2 QSLs. 

Beyond mean field theory, apart from the gap- 
less/gapped spinon excitations, there exist additional 
low-energy excitations, depending on the gauge struc- 
ture. For example, in gapless U(l) QSL in two spatial 
dimensions, there is an emergent gapless photon. For Z2 
spin liquids, the "magnetic flux" associated with the Z2 



QSL Label 


Singlet Ansatz 10 


"Singlet + Triplet" Ansatz 


C7(l)[0,0] 


Fermi surface (F.S.) 


F.S. is altered 
but not gapped out 


tf(l)[0,7T] 


Dirac point (D.Pt.) 


D. Pt. is gapped out (*) 


C7(1)[7T,7T] 




The Flat band 




Flat bands 


acquires dispersion 


E7(1)[tT,0] 


& D.Pt. 


while the D. Pt. 
remains intact. 


Z 2 [0, 0]A 


F.S. gapped out to 


Band touching points 




Band touching points 


remain 


Z 2 [0,tt]/3 


F. S. becomes 


Gap is 




fully gapped 


altered 


Z 2 [0,0]B 


F.S. shrinks compared 


F.S. gapped out to 




to parent U(l) 


band touching points 


Z 2 [0,7i]a 


D.Pt. changes to 


Band touching points 




band touching points 


remain 


Z 2 [0,0]D 


F.S. shifted compared 


F.S. gapped out to 




to parent U(l) 


band touching points 


Z 2 [0,7T] 7 


D.Pt. changes to 


Band touching points 




band touching points 


remain 


Z 2 [tT,7t]B 


Negligible change 


Bands gapped to D.Pt. 




compared to 


and band touching points 




parent U(l) 




Z 2 [tt,0]B 


Flat Bands are gapped 


F.S. gapped to 




to form a F. S. 


band touching point 



TABLE I. (* This state is unstable to confinement). Struc- 
ture of low-energy spinon excitations in the U(l) and Z 2 under 
consideration, in the mean-field parameter regime described 
in Sec. IIIB. These behaviours persist within the neighbour- 
hood of the respective parameter sets; it gives a general under- 
standing of these states. The singlet case shows the evolution 
from £7(1) to Z 2 spin liquids upon addition of the bolded terms 
in Table II. The "singlet + triplet" case shows the evolution 
of all these states upon addition of nearest-neighbour triplet 
terms. For Z 2 spin liquids, F.S. indicates a Fermi surface of 
Fermionic quasi-particles, obtained upon diagonalization of 
the mean- field Hamiltonian 1~Lq. The low energy excitations 
lie about the spinon Fermi-level which is set by the half filling 
constraint on the spinons (Eq. 6). 



gauge field leads to gapped topological excitations called 
visons. 63 In this work, we shall be mostly concerned with 
the question of the nature of the ground state and spinful 
excitations, such as spinons and spin-1 excitations. 

Having identified the spin liquids, we turn to their dy- 
namical spin-structure factor, which is just the Fourier 
transform of spin-spin correlation function. The dynam- 
ical spin-structure factor has been recently measured in 
inelastic neutron scattering experiments on a single crys- 
tal sample of Herbertsmithite. It shows diffuse scattering 
with no signature of the spin gap down to at least 1.5 
meV. 9 

In our computations, almost all of the above spin-liquid 
ansatze show patterns of diffuse intensity in the dynam- 
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ical structure factor throughout the Brillouin zone, the 
overall shape of which is largely influenced by the parent 
U(l) spin liquid. The structure factor of the U(l) states 
have characteristic flat and dispersive features within the 
diffuse continuum. In the Z 2 spin liquids, the dispersive 
features are broadened and become significantly more 
complex. A central result is that, for some Z^ spin liq- 
uids, these dispersive features are only slightly stronger 
than the diffuse background, so that the dynamical struc- 
ture factor looks mostly diffuse. This is particularly true 
for the states labelled Z2[0,7r]/3 and Z2[0,7r]a (the dy- 
namical structure factor is shown later in Figs. 10 (e) and 
(f), respectively). This may be consistent with present 
neutron scattering results where similar diffuse scattering 
is seen over a large energy window in large parts of the 
Brillouin zone. 9 It is found that the spin-rotation symme- 
try breaking terms primarily introduce non-zero intensity 
at the centre of the Brillouin zone. They also split the 
spinon bands, so that diffuse scattering is found more 
evenly throughout all of the Brillouin zone in the pres- 
ence of DM and Ising interactions. The inelastic neutron 
scattering indeed shows some intensity at the centre of 
the Brillouin zone, with a broad maximum near 0.4J. 

To better understand the effect of triplet terms on 
the ansatze we consider the ESR absorption spectra 
that show a non-trivial response solely from spin-rotation 
breaking terms. The ESR signal shows the largest 
absorption intensity at a (5-function peak at the Zee- 
man field energy scale. Triplet terms create additional 
"satellite" absorption peaks, offset by their energy scale. 
The number, position, and broadening of these ex- 
tra absorption peaks differentiate between the different 
ansatze. Particularly for the aforementioned Z2[0,7r]/3 
and Z2[0,7r]a states, the absorption spectrum provides 
a distinct qualitative features (as in Figs. 14 (e) and (f)) 
that can be verified experimentally. 

The rest of this paper is organized as follows. We 
begin in Sec. II by discussing the spin model used to 
capture the behaviour of Herbertsmithite. Following the 
introduction of the spin Hamiltonian, in Sec. Ill we de- 
scribe the fermionic slave-particle construction of mean- 
field states. In Sec. IV we describe the projective sym- 
metry group analysis of the Z2 ansatze, describing their 
respective U(l) parent QSL states as well as the effect 
of the spin-rotation symmetry breaking terms. The spin 
correlations of these states are characterized in in Sec. V. 
Finally, in Sec. VE we discuss the implications of these 
correlations in light of the recent neutron scattering re- 
sults. We summarize our results in Sec. VI. The details 
of different calculations are given in various appendices. 



II. SPIN-1/2 MODEL ON THE KAGOME 
LATTICE 

The starting point for our analysis is the spin-1/2 
model with antiferromagnetic nearest-neighbour Heisen- 




FIG. 1. (Color online). Bond labels and directions for nearest- 
neighbour (dark blue, solid) and next-nearest-neighbour (light 
brown, dashed) bonds on the Kagome lattice. The nearest- 
neighbour directions also indicate the directions of the 
Dzyaloshinsky-Moriya terms. The bonds labelled 1 are the 
ones upon which mean-field parameters are defined as in Ta- 
ble II; parameters on other bonds are determined from these 
through the application of symmetry operations. 




FIG. 2. (Color online). Symmetry transformations of the 
Kagome lattice, including translations T\ and T 2 , a rotation 
Cq and reflection cr. 



berg interactions on the Kagome lattice, 

^nn-khaf = J ^ &i ' Sj. (1) 

where J ~ 17 meV ~ 200 K denotes strong antiferro- 
magnetic interactions, and forms the largest energy scale 
in Herbertsmithite. 14,15 Experiments also indicate that 
there may be three kinds of interactions acting as per- 
turbations to the above Hamiltonian: 

a. Dzyaloshinsky-Moriya (DM) interaction: 

Hvu = Dj2Vij' s i x S r (2) 

This has been estimated to have energy scale of |D| ~ 
15 K (perpendicular to the Kagome planes) from ESR 
experiments 13 and fitting of the anisotropy of the mag- 
netic susceptibility data. 8 In Eq. (2), the orientation of 
T>ij = z is out of the plane when the bonds are counted in 
a counter-clock wise way around each triangle (Fig. 1). 
The in-plane components of the DM interactions (~ 2 
K) appear to be negligible in comparison with the out- 
of-plane component, and we will neglect it in our calcu- 
lations. 

b. XX Z anisotropy: The susceptibility anisotropy 
also indicates a sizable easy-axis spin-spin interaction 
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along an axis perpendicular to the Kagome plane, 

Rising = A ^2 SiSji ( 3 ) 

m 

where A - J/10 - 20K. 

c. Next-nearest-neighbour coupling: As discussed in 
the introduction, even small J 2 can drastically affect the 
stability of spin liquids of the HKAF, as found in re- 
cent numerical calculations. Accordingly, we consider 
adding an isotropic next-nearest-neighbour antiferromag- 
netic Heisenberg interaction J2. 

^NNN = J2 ^ &i ' Sj. (4) 

The complete Hamiltonian is then obtained by adding 
Eqns. (1) - (4). All these terms have full symmetries 
of the Kagome lattice, shown in Fig. 2, and enumerated 
below. These are: 

• Translations Ti and T2 in the Kagome plane. 

• A sixfold rotation Cq around the centre of a 
hexagon. 

• A reflection a through a hexagon. 

With this Hamiltonian, we now construct possible Z2 
spin liquid states with fermionic spinons, both gapped 
and gapless, and then contrast their properties in con- 
text of experiments on Herbert smithite. We shall limit 
ourselves to Z2 spin liquid ansatze that can be realized 
with just first and second-neighbour terms (both hopping 
and pairing). Since mean-field theory only gives an order 
of magnitude estimate of the microscopic parameters, we 
shall investigate the existence of the above spin liquids in 
somewhat a broad range of the parameters, especially J2, 
whose effect is known to be quite sensitive in stabilizing 
Z2 spin liquids. 7 



III. SLAVE- FERMI ON CONSTRUCTION OF 
QUANTUM SPIN LIQUID STATES 

A. Spinon Representation 

In the slave- fermion formalism, as noted earlier, we 
represent a spin- 1/2 operator, at a site z, in terms of a 
fermion (spinon) blinear: 

where a,/3 G {t?i} are the two flavours of spinons. 54 57 
Here, cr a (a = 1, 2, 3) are the Pauli matrices. The Hilbert 
space of the spinons is twice as large as the original S = 
1/2 Hilbert space. The physical spin wave- functions are 
obtained by projecting the spinon wave-function in the 



subspace of single spinon per site. This is equivalent of 
having the following constraints for every lattice site: 

fiafia ~ 1 

fiafi{3ea{3 = = 0> (6) 

where 612 = —621 = 1 is the 2x2 completely antisym- 
metric tensor. 



B. Mean-Field Decoupling 

With the above spinon representation, we can re-write 
the bilinear spin Hamiltonian of Sec. II as a quartic 
fermion Hamiltonian in terms of the spinon operators. 
We write the Heisenberg, DM and Ising anisotropic ex- 
change interactions on the nearest neighbour bonds as 
the following: 

+ ^(-3'?l i %+Y] i -Y ij ), (7) 

z-S i xS j = - Etjxij) + ffijY* ~ YgfUj), 

(8) 

SIS', = ^(-xtjXij - + + E\}El) 

+ " + + Y u Y ^ ( 9 ) 

We have defined 

Xij = fiafiapfjP, = fi a T%pfj/3 

Vij = fia(iT 2 ) a pfj(3, Y°j = j \ a (iT* V a ) a/3 j * j(3 m (10) 

(a = x, y : z) the singlet and triplet hopping (particle- 
hole) channels and singlet and triplet pairing (particle- 
particle) channels, respectively. Note that unlike the 
usual slave fermion decomposition, where only the singlet 
decoupling channels are used, we have introduced both 
singlet and triplet decoupling channels. As discussed ear- 
lier, since spin-rotation symmetry is already broken in the 
Hamiltonian, there is no a-priori reason to choose only 
the singlet or triplet channel. Indeed, in the following we 
shall see that the triplet channel plays an important role 
in our calculations. 

These quartic operators are then decomposed by intro- 
ducing eight auxiliary fields corresponding to E and 
Y for all the particle- hole and particle-particle channels. 
The saddle point approximation in terms of these aux- 
iliary fields then gives the mean-field quadratic fermion 
Hamiltonian, 1-Lq, characterized by the non-zero decou- 
pling channels. 

In the singlet case, %q is conveniently written using 
a two-component Nambu basis which most clearly 
reveals the underlying gauge redundancy of Hq. The 
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Hamiltonian is given by 28 ' 29 ' 56 ' 58 

ij 

where the Uij matrix satisfies Uji = Ujj, and 



(ii) 



TrrSing 




(12) 



The exact constraint of single spinon per site (Eq (6)) 
is then relaxed to an average constraint and imposed by 
Lagrange multipliers, given by 28 ' 29,56 ' 58 



(13) 



where \i a (a = 0, 1, 2, 3) are the Lagrange multipliers. We 
note that the Nambu (z^) basis components involve op- 
erators of the same spin type. This form clearly demon- 
strates the explicit SU (2) gauge redundancy 28 ' 29,56 ' 58 



VWi e SU(2). 



(14) 



The introduction of triplet terms necessitates, in prin- 
ciple, a four-component basis. 59 ' 62 However, since in the 
following we will use only the z-components of the triplet 
terms, we may use the same basis Vi as the singlet case, 
writing terms with the matrix form L^ trip , given as 



U 



trip 



El 



-Y z E z - 

13 ij j 



(15) 



The quadratic form for the triplet part is now given by 
v}Uij ip Vj. Importantly, this form has the same gauge re- 
dundancy as the singlet case, so the projective symmetry 
group classification will be similar, as to be discussed in 
Sec. IV. These terms combine in the quadratic Hamilto- 
nian to give 



n Q = j2»hu:p+u^)v j+ h.c. 



(16) 



As pointed out earlier, with the introduction of the DM 
and Ising anisotropy terms on the NN bonds, the singlet 
and triplet channels mix and the saddle point is char- 
acterized by a combination of the two kinds of channels. 
For example, in the hopping sector, we can re-diagonalize 
the interactions to write (for NN bonds): 



J^i ' Sj 



W-X-X- 



- W+X+X+ 



f as? s; ■+ 

—E x] E x 



Dz ■ S,- x S, 



J- 



—Ey^E y , 



(17) 



where we withhold the z, j subscripts on all operators for 
the rest of this section, and 



J + A , 1 



o 4 

VJ 2 +Z) 2 T J , . D * 
X± ~, X T i-z-E , 

|^ ± | 2 = 2 V / J 2 T^ 2 (VJ 2 + =p j) . (18) 

For the pairing channels, the same is true with rj in place 
of and Y in place of E. 

Both x- and ^-components of the triplet terms lead 
to unstable mean-field states, so we will only consider 
ansatze with singlet (x, rj) and z-component triplet (E z , 
Y z ) terms. We can write the mean-field Hamiltonian in 
the basis described by (12) and (15). On a NN bond, we 
have 

«g N = w_ (x*_X- + V-V-) + ^+ (x* + x+ + V+V+) 

+h.c. - oj-\x-\ 2 + \v-\ 2 - uj+Ix+I 2 + \v+\ 2 - (19) 

On a NNN bond (since there is only antiferromagnetic 
Heisenberg term, and no D or A) we decouple only in 
the singlet channels: 



(20) 



/NNN 



T-Uy^^ can now 



The quadratic Hamiltonian Hij = 

be solved to obtain the spinon wave- function, |^ /spinon . 

To gain insight on the actual values of different mean- 
field parameters that characterize the different QSL 
ansatz, we perform self-consistent mean-field calcula- 
tions, the details of which are outlined in Appendix C. 
We note that in the self-consistent calculation, we find 
that the next-nearest-neighbour exchange, J2, necessary 
to stabilize six out of the eight Z2 spin liquids in Table II 
is: J2 > 0.4J. This large value of J2, required to stabilize 
non-zero next nearest neighbour hopping and pairing in- 
teractions for the spinons, is certainly an artifact of our 
mean field treatment. 

These mean-field results leading to stable solutions 
act as a guide for choosing a representative mean-field 
parameter set for calculation of the dynamical spin- 
structure factor. The general qualitative features of the 
structure factor so obtained are expected to be less sensi- 
tive to the details of the exact parameter values. In any 
case it is likely that the parameter values obtained from a 
mean-field theory are renormalized by quantum fluctua- 
tions. Hence, we take mean-field parameters of singlet Z2 
terms wherein pairing and hopping terms are of similar 
magnitudes, to demonstrate the qualitative features of 
these Z2 states in comparison to the U(l) ones. We will 
also take next-nearest-neighbour terms on the order of 
0.4 times the magnitude of the nearest-neighbour terms, 
and a triplet-singlet ratio of 0.1 ~ D/J when considering 
a spin-rotation symmetry-breaking ansatz. 



7 



IV. PROJECTIVE SYMMETRY GROUP AND 
MEAN-FIELD ANSATZE 

The spinon wave- function \^) sp i n on so obtained is not a 
valid spin wave-function, since it only satisfies constraint 
Eq. (6) on average. Hence, it must be projected into a 
single-spinon-per-site subspace to obtain the spin wave- 
function through numerical Gutzwiller projection. 28 

Alternatively, it is possible to consider the low energy 
fluctuations about the mean-field solution. The gauge 
redundancy of the spinon description (Eq. (14)) imme- 
diately suggests that the low energy effective theory con- 
tains spinons minimally coupled to a lattice gauge field. 
The structure of the gauge group can be determined from 
the structure of the mean field QSL ansatz for T~Lq 28,62 
in the following way. Each mean-field ansatz is char- 
acterized by Uij G SU{2) on the bonds (Eq. (12) and 
(15)). Taking the product of such UijS over various closed 
loops (C) starting from a specific base site (i), we obtain 
the "flux" of the Uij link fields which has the following 
form: 29 

WcW= II u *= E A c^ Ta ( 21 ) 

(jk)eC a=0,x,y,z 

where Aq(i) are numbers specific to the loop and r a are 
the identity (a = 0) and Pauli matrices (a = 1,2,3) in 
the gauge space. For U(l) QSLs, A a (i) for different loops 
are proportional to each other for the same a and at least 
one of A a (i) is non-zero for a = 1, 2, 3 (say 3). Then it is 
possible to write 29 

C/(1)QSL : W c (i) oc r 3 VC (22) 

This means that there exists a gauge in which all Uij oc 
t 3 and hence Hq of the form in Eq. (11) is invariant 
under all gauge transformation: Vi —> e 10 ^ 3 ^ where, 0i G 
[0, 2ir]. Thus the gauge group is U(l) and accordingly we 
have a U(l) QSL. 29 

For Z 2 spin liquids, two or more loops based at the 
same site have Aq(i) that are not proportional to each 
other and a similar analysis as above shows that the 
gauge transformation of the following form is only pos- 
sible: Vi — ^ Q^i, where Q = ±1. 29 This is a Z 2 gauge 
group describing Z 2 spin liquids. The above argument 
can be extended to the case of the most general triplet 
ansatz with a four- component spinon representation. 62 

The gauge group so obtained is called the invariant 
gauge group, (IGG) because the mean field ansatz re- 
mains invariant under such gauge transformations. Due 
to the above gauge redundancy, the lattice symmetries 
(and time reversal) under which the Hamiltonian is in- 
variant act project ively on the mean-field ansatze. This 
means that usual physical symmetry transformations on 
the ansatz transforms it to its gauge equivalent form; 
hence an added set of gauge transformation is required 
following the actual symmetry transformation to bring 
the ansatz back to itself. The PSG characterizes, for a 
given set of symmetry transformations SG, the possible 



distinct sets of gauge transformations {Gs} for each op- 
eration S G SG, where 

(GsS)Uij(GsSy = Uij. (23) 

A particular set {Gs} characterizes a corresponding 7-Lq 
and hence the mean field ground state. Therefore, ele- 
ments of the IGG are the gauge transformations Gj as- 
sociated with the identity operation I of the symmetry 
group SG. In other words, the symmetry group SG = 
PSG/IGG. One can choose a gauge where the elements 
of the IGG are site-independent, namely G/(i) = ±7 for 
Z 2 transformations, Gj(i) = e 10 ^ (0 G [0,2tt)) for U(l) 
transformations, and G/(i) = e %dn '^ (n G S 2 , 6 G [0, 2tt)) 
for SU(2) transformations, characterizing Z 2 , U(l) or 
SU(2) spin liquids, respectively. 29 

Determination of the allowed sets {Gs}, and thereafter 
the ansatze Hq, comes from the group multiplication ta- 
ble of SG. For every group multiplication rule written 
in the form ABC -1 = 7, a similar constraint is put 
on the PSG, namely (G A A)(G B B)(G C C)- 1 G IGG, a 
pure gauge transformation. These constraints determine 
the allowed PSG {GsS} for the symmetry group of the 
Kagome lattice. The singlet Z 2 spin liquids have been de- 
termined by Lu et. a/., 10 and the same PSGs apply in the 
triplet case, since we use the same basis (12). However, 
the resulting ansatze differ between these cases, primar- 
ily in the existence and structure of the triplet terms. We 
will leave these details to Appendix. B, and summarize 
the results here. 



A. Singlet Mean-Field Ansatze 

In this section, we characterize the singlet spin liq- 
uid ansatze on the Kagome lattice that are of interest to 
us. 10 These include the eight Z 2 spin liquids that occur 
in the vicinity of one of four different parent U(l) spin 
liquids and can be stabilised by tuning spinon hopping 
and pairing terms up to the second nearest neighbour 
sites. We shall begin by briefly describing the spinon ex- 
citation spectra of the four parent U(l) spin liquids at 
the nearest-neighbour level. The Z 2 states are obtained 
upon appropriate introduction of additional hopping and 
pairing terms that break the IGG to Z 2 , as are shown in 
Table II. 

The nomenclature of the four U(l) QSLs have al- 
ready been introduced in Section I A. These are iden- 
tified as 17(1) [0,0], E/(1)[0,tt], 17(1) [tt, tt], and C7(1)[tt, 0] 
phases. In Appendix A we show the signs of the hop- 
ping parameters used (in the gauge that we have chosen) . 
Here we summarize the features of their spinon spectrum 
in the above gauge. We note once again that though 
gauge dependent, the spinon spectra provides valuable 
clues to the nature of the experimentally measurable 
and gauge invariant spin-spin correlation functions. We 
plot the dispersion along the high symmetry direction 
T^M^K^Toi the original Brillouin zone (Fig. 3). 
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FIG. 3. The cut along the high symmetry direction of the 
primitive and extended first Brillouin zone is shown respec- 
tively. The boundary of the primitive (extended) Brillouin 
zone is shown in light (dark) gray. While the spinon disper- 
sion is periodic within the primitive Brillouin zone, the spin- 
structure factor is periodic in the extended Brillouin zone (see 
text for details). 
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FIG. 4. (Color online). Spinon dispersion for the (a) 
C/(1)[0,0], (b) C/(1)[0,tt], (c) I7(l)[7r,7r], and (d) U(l)[w,0] 
singlet spin liquid ansatze. These dispersion go from the cen- 
tre of the Brillouin zone of the Kagome lattice, T, to the edge 
M, to the corner K, and back to the center T. The Fermi 
level is at zero energy. 



J7(l)[0,0]; The uniform [7(1) [0,0] state has a Fermi 
surface of spinons. Since this state has uniform hop- 
ping, its dispersion is simply the band structure for the 
Kagome lattice, and is seen in Fig. 4 (a). It has a fiat 
band at the maximum of the dispersion. 

J7(l)[0,7r]; This ansatz breaks translational symme- 
try of the Kagome lattice, and has Dirac points in its dis- 
persion, as seen in Fig. 4 (b). It has a doubly-degenerate 
flat band at the maximum of the dispersion, similar to 
the U(l) [0,0] case. 

U{1) [it, 7r] ; This state has a Dirac point at the centre 
of the Brillouin zone, but also has double-degenerate flat 
bands at the Fermi level, as seen in Fig. 4 (c). 

U(i)[ir, 0]: We find that there are at least two dis- 
tinct U(l) states with a [7r, 0] hopping flux pattern. Only 
one of them can stabilize a Z 2 spin liquid {Z^, 0]B) 
upon addition of the appropriate next-nearest-neighbour 
hopping and pairing terms; therefore, this is the state 
that we will consider. Like the J7(l)[7r, 7r] state, there 
are also flat bands at zero energy, as seen in Fig. 4 (d). 



SL Label 


Unn 


Unnn 


Z 2 [0, 0]A 


2 .3 

T ,T 


T 3 


Z 2 [0,tt]/3 


2 .3 
T ,T 


T 3 


Z 2 [0,0]B 


.3 

T ,T 


r 3 


Z 2 [0,7i]a 


2 3 

T ,r 


r 3 


Z 2 [0,0]D 


T 3 


T ,T 


Z 2 [0,tt]j 


r 3 


r 2 ,r 3 


Z 2 [k^]B 


r 2 


T 3 


Z 2 [tt,0]B 


r 2 


T 3 



TABLE II. The general structure of the singlet QSL ansatz 
for the Z 2 spin liquids. The particular bonds chosen for Unn 
and Unnn are shown in Fig. 1. The denoted terms in Unn 
and Unnn show the form of the pairing and hopping terms 
required to break the U(l) parent spin liquid's gauge group 
down to Z 2 . We note that we have chosen a gauge in which 
t 2 terms for Z2[0,0]A and Z 2 [0,tt]/3 ansatze are zero. 10 



However, there is a gap to the subsequent single-spinon 
excitations throughout the entire Brillouin zone. 

On adding appropriate spinon hopping and pairing 
terms up to second nearest neighbour, the IGG of the 
above spin liquids is broken down from U(l) to Z 2 . 10 
There are eight such Z2 spin liquids as given in Table II. 
The addition of these terms changes the structure of the 
low-energy quasiparticle excitations, with parameters as 
discussed in Sec. IIIB. In some cases, a gap is opened 
(Z 2 [0, 7r]/3). In others, the line degeneracy of zero-energy 
excitations from the U(l) state's Fermi surface remains 
(^ 2 [0, 0]B and Z 2 [0, 0]D). For some other states, only 
band touching points remain (Z 2 [0,0]A, Z 2 [0,7r]a and 
Z 2 [0, The flat bands at zero energy of the J7(l)[7r, tt] 
and J7(l)[7r, 0] states either remain (Z<i\k, ir\B) or are 
lifted to a line degeneracy {Z<£(k, 0]B). Table I summa- 
rizes the structure of the low-energy quasiparticles for the 
above U(l) and Z 2 spin liquids. The full form of {Gs(i)} 
is given in Appendix B. 



B. "Singlet + Triplet" Mean-Field Ansatze 

All of the above spin liquids on the Kagome lattice are 
also realized even in absence of spin-rotational symmetry, 
if time-reversal and the symmetries of the Kagome lattice 
are preserved. The allowed Uij now has both real singlet 
terms and imaginary z-component triplet terms for the 
nearest neighbour hopping and pairing, and real singlet 
terms for next-nearest-neighbour as well. 

Again, we first consider the effect of the triplet terms 
on the four U(l) spin liquids. In all cases, the spec- 
tra show additional features in their dispersions, com- 
ing from the triplet terms that removes the flat bands 
previously present throughout the Brillouin zone. The 
C/(1)[0, 0] state has an altered Fermi surface, which does 
not split the fa and f± spinons (Fig. 5 (a)). The Dirac 
cone in the J7(l)[0,7r] state, (Fig. 4 (b)) however, is 
gapped out (Fig. 5 (b)). Further, the bands are also 
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FIG. 5. (Color online). Spinon dispersion for the (a) 
17(1) [0,0], (b) f7(l)[0,7r], (c) C7(l)[7r,7r], and (d) *7(1)[tt,0] 
spin liquid ansatze with "singlet + triplet" channels for NN 
bonds. These dispersion go from the centre of the Brillouin 
zone of the Kagome lattice, T, to the edge M, to the corner 
K, and back to the center F. The fermi level always lies at 
the zero of the energy scale. We note that since the "singlet 
+ triplet" C7(l)[0,7r] spin liquid state is gapped, it is unstable 
to confinement due to instanton tunnelling events. 64 



spin split throughout the Brillouin zone. We note that 
such a two-dimensional gapped [7(1) spin liquid is unsta- 
ble toward a confinement transition. 39,64 So we infer that 
triplet decoupling may indirectly render the [7(l)[0,7r] 
state unstable. The U(1)[tt, it] and [7(l)[7r,0] states both 
have their zero-energy flat bands mostly gapped out, 
leaving a filled band with a line degeneracy at zero en- 
ergy. The Dirac cone in the [7(l)[7r,7r] state remains in- 
tact (5 (c)), while newer Dirac nodes are created for the 
[7(l)[7r,0] state on including the triplet terms (Fig. 5 
(d)). 

Allowing for small triplet terms in the singlet Z^ 
ansatze also changes the low-energy structure of the 
quasiparticles. The Z2[0,7r]/3 state has its gap altered. 
For states with line degeneracies, we find that only 
band touching points remain (^[0, 0]B, Z2[0,0]D and 
Z2[tt,0]B). Also, the band touching points are found to 
be robust to spin-rotation symmetry breaking perturba- 
tions (Z 2 [0,0]A, Z 2 [0,7r]a) and Z 2 [0,tt]j). In the state 
^[tt? k]B, the flat bands at zero energy are lifted to band 
touching points, while the Dirac node persists. These 
changes are summarized in Table I. 

With this, we have completed the description of the 
candidate QSL states. We shall now calculate the dy- 
namical spin-structure factor and study their broad fea- 
tures that can help us identifying the nature of the spin 
liquid, possibly realized in Herbertsmithite, from inelas- 
tic neutron scattering and ESR experiments. 



V. SPIN CORRELATIONS IN U(l) AND Z 2 SPIN 
LIQUIDS 

A. Structure Factor and Experimental Probes 

To determine the spin correlations, we will calculate 
the dynamical structure factor; the matrix is given by 



-c 



dt 



Auot 



5^ e *-('«-'i)( S ?(t)Sj(0)), 

(24) 



where a,/3 G {1,2,3}. It characterizes spin-1 magnetic 
excitations of energy u and wavevector q. Up to a mag- 
netic form factor, S af3 (^u) can be measured directly 
by inelastic neutron scattering. 9 We note that Eq. (24) 
is not periodic in the reciprocal lattice vectors of the 
Kagome lattice. Since Ti — rj can be half of the primitive 
lattice vectors a or b in Fig. 2, we must extend the Bril- 
louin zone to double its size. We will plot the dynamical 
structure factor along the path T — > M' — >> K' — ^ T in 
the extended Brillouin zone, where T is the centre, M' 
is the midpoint of an edge, and K' is the corner of the 
edge of the extended Brillouin zone (EBZ) (Fig. 3). We 
found that the equal-time (cj-integrated) structure factor 
shows little qualitative differences among our ansatze, so 
we concentrate on the unresolved features for different 
QSL states. Further, for the singlet ansatze, since the 
quadratic Hamiltonian Hq commutes with the total spin 
operator and the ground state is an eigenstate of the to- 
tal spin with eigenvalue 0. Thus, the dynamical structure 
factor of singlet ansatze is zero at q = 0. 

The dynamical structure factor can also be probed by 
ESR, under an external Zeeman field Hz J2i &i ' Ab- 
sorption of a transverse microwave field (along axis a) of 
frequency uo probes the long- wavelength q w dynamical 
structure factor S aa (0, u). 49 The absorption intensity is 
given by 



I(w) = 



H^ 



S aa (0,co), 



(25) 



where the field is applied in the a direction, with ampli- 
tude H m and frequency uj. 50 The SU( 2) -invariant terms 
affect the intensity trivially; that is, I(uu) oc S(u — Hz), 
where Hz is the Zeeman field strength. Hence, ESR line- 
shape is sensitive to spin anisotropy, and so the triplet 
terms. Thus, ESR experiments can reveal important in- 
formation about the effect of the spin-rotation symmetry 
breaking perturbations. 



B. Dynamical Structure Factor of Singlet U(l) Spin 
Liquid Ansatze 

We will begin by considering the singlet U(l) spinon 
ansatze. Our choice of color scale is made to accentu- 
ate the low-intensity scattering compared to the zero- 
intensity background. 
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FIG. 6. (Color online). Dynamical structure factor of the 
NN singlet £7(1) [0, 0] spin liquid along the high symmetry di- 
rection of the extended Brillouin zone. The characteristic low 
energy dome-like structure results from the spinon excitations 
near the fermi surface (Fig. 4 (a)). 



FIG. 7. (Color online). Dynamical structure factor of the 
nearest-neighbour singlet £7(l)[0,7r] spin liquid along high 
symmetry directions of EBZ. Compared to [7(1) [0,0] state, 
the low energy continuum near T and M' points are replaced 
with cones of scattering which is a consequence of the Dirac 
node at the fermi level for the spinons (Fig. 4 (b)). 
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FIG. 8. (Color online). Dynamical structure factor of the 
nearest-neighbour singlet U(1)[tt, tt] spin liquid along the high 
symmetry direction of the EBZ. The extensive continuum of 
low energy scattering is a contribution of the flat bands at the 
fermi- level (Fig. 4 (c)). 



FIG. 9. (Color online). Dynamical structure factor of the 
nearest-neighbour singlet £7(l)[7r,0] spin liquid along high 
symmetry directions of the EBZ. There is finite scattering 
intensity exactly at zero energy because of the flat bands in 
the spinon dispersion (Fig. 4 (d)), which is followed by a lack 
of scattering due a to gap in the spinon spectrum. 



Fig. 6 shows the dynamical spin-structure factor for 
the £7(1) [0,0] state. The low-energy domes of scattering 
anchored around the T and M' points are contributed 
by the excitations near the Fermi surface. Above these 
domes, we see additional flat scattering coming from con- 
tribution of the the higher energy, fiat bands (refer to Fig. 
4(a)). 

In Fig. 7 we plot the dynamical spin-structure factor 
for the J7(l)[0,7r] state. In contrast with the 17(1) [0,0] 
case, low-energy cones of scattering are seen at the T, 
M' , and intermediate points in the extended Brillouin 
zone, a consequence of the linearly-dispersing low-energy 
single-spinon excitations of the Dirac spin liquid. The 
flat, intense band of scattering seen at the highest ener- 
gies are manifestations of the flat spinon bands (Fig. 4 
(b)). 

Fig. 8 shows the structure factor for the C7(l)[7r,7r] 



state. In this case, we can see intensity all the way down 
to zero energy across the extended Brillouin zone. Within 
this diffuse scattering, there are two major variations of 
intensity. The first is dispersive, rising from the T point, 
due to the Dirac point in the spinon dispersion. The 
second is the flat band of intensity in the middle of the 
energy range, again due to zero-energy flat bands of the 
dispersion. 

Finally, Fig. 9 shows the structure factor for the 
£7(1) [7r, 0] state. While the flat bands in the dispersion 
contributes to the scattering exactly at zero-energy, the 
gap to subsequent excitations is seen in the absence of 
scattering. The highest intensity is seen at a particular 
point of scattering at M' . 

We see that the low-energy spin correlations are an 
effective way to distinguish between these C7 (1) spin liq- 
uids. Other dispersive scattering features characteristic 
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to the states also show up at comparatively higher ener- 
gies. 

Next, we will look at the effect of the additional pairing 
terms of the Z^ spin liquid ansatze on the dynamical 
structure factor, particularly as a means to distinguish 
between Z2 states with the same parent U(l) spin liquid, 
in spite of their similarities. 



C. Dynamical Structure Factor of Singlet Z2 Spin 
Liquid Ansatze 

We see that the dynamical structure factor of Z2 states 
are similar in overall shape to that of their parent U(l) 
states. To illustrate this point, we will consider the 
Z 2 [0,0]A, B and D states, in Fig. 10 (a,b,c). The 
^[0, 0]D state has a structure factor that is most similar 
to the parent £7(1) [0,0] case (Fig. 6), with zero-energy 
scattering along much of the M' — K' line. Above that, 
the scattering is diffuse, with dispersive features that are 
only slightly stronger than the diffuse background. For 
the ^[0,0] A state, the low-energy domes of scattering 
persist over small energies, and the flat features have 
become much more dispersive. There is a large inten- 
sity at the M' point at higher energies. In the Z2[0,0]£? 
case, there is low-energy intensity only near the T and 
M' points. Many dispersive features are seen, particu- 
larly near the K' point. 

The Z 2 [K 1 'n)B state has a dynamical structure factor 
(Fig. 10 (d)) which is strikingly similar to its parent 
J7 (1) [7r, 7r] state (Fig. 8). Only the intensity of the dis- 
persive bands increases, and scattering close to T is found 
up to all energies, but there are no obvious qualitative 
differences between the two. 

We also see similar resemblance in case of the the 
Z2[0,7r] QSL states with their parent i7(l)[0,7r] state 
(Fig. 10 (e,f,g)). The Z2[0,7r] states retain flat inten- 
sity at the highest energy, and relatively constant diffuse 
scattering across their energy range from their parent 
state. The Z2[0,7r]/3 state has a fairly significant spin- 
gap, where the low-energy cones of scattering of the par- 
ent U(l) state (Fig. 7) have been rounded off. This spin- 
gap is found to be rather momentum independent. Also 
the flat features of the £7(1) [0, n] state are found to be 
broadened, and the scattering is fairly diffuse, although 
some dispersive intensity can be seen coming from the 
T point, due to the remnant of the Dirac cone. At low 
energies, the intensity for the Z2[0,7r]a state has many 
broad domes of scattering, and zero-energy excitations 
exist almost throughout. At larger energies, the intensity 
is mostly diffuse, where the C/(l)[0,7r] state's high-energy 
intensity has been split into two neighbouring broad, dif- 
fuse, yet still almost flat bands. The low energy features 
of the Z2[0,7r]7 state is similar to the U(l) case. Above 
this, the intensity is also separated into two broad, diffuse 
bands, more prominently than the Z2[0,7r]a state. 

The Z2[7r,0]B state (Fig. 10 (h)) shows many sepa- 
rated broad and diffuse bands, including the band near 



zero energy, which is the remnant of the strictly zero 
energy scattering coming from the flat bands of the 
f7 (1) [7r, 0] state (see Fig. 9). Furthermore, the intensity 
is relatively constant across the cut, only diminishing as 
the zone centre (r point) is approached. 

With the exception of the [-zr, tt] states, the low-energy 
intensity and dispersive features allow the twelve spin- 
liquid ansatze in Table II to be distinguished from each 
other qualitatively. In the next subsection, we will see 
the effect of adding triplet terms to the QSL ansatze to 
the dynamical structure factor. 



D. Dynamical Structure Factor of "Singlet + 
Triplet" U(l) and Z2 Spin Liquid Ansatze 

The introduction of the triplet terms E\- and 
breaks the spin-rotational symmetry of the spin-polarized 
dynamical structure factor S a ^(q_, cj), where S xx = 
gyy _^ gzz _ while the singlet structure factors showed 
a vanishing intensity approaching q = 0, these triplet 
terms generate a non-zero intensity in S xx and S yy . Also, 
the scattering intensity at low energies changes, due to 
changes in the spinon dispersion (summarized in Table. 
I). Since the spinon bands are, in general, split by triplet 
terms, the scattering becomes increasingly diffuse. 

We will showcase these changes by calculating the trace 
of the dynamical structure factor matrix, i.e., (S xx + 
S yy + S zz )/3, and then plotting the difference from the 
singlet structure factor. In Appendix. D, we show the 
actual structure factor for all these "singlet + triplet" 
states. 

The "singlet + triplet" J7(l)[0, 0] structure factor looks 
very similar to the singlet case, particularly in the low- 
energy scattering. However, the diffuse scattering at 
moderate energies is no longer flat. While this is also true 
for the "singlet + triplet" C/(1)[0, 7r] structure factor, it 
also has two qualitative differences. The first is that the 
cones of scattering at low energy are gapped out; this 
can be seen more clearly in the full "singlet + triplet" 
structure factor in Appendix D, Fig. 19. The second is 
that the flat band at the highest energy is split into two. 
For the J7(l)[0,7r] state, though we plot the structure 
factor for completeness, we note that in this state the 
spinon spectrum is gapped and hence the state is unsta- 
ble to confinement transition due to instant on events. 64 
The "singlet + triplet" J7(l)[7r, 7r] structure factor is very 
similar to the singlet case, where the low-energy scatter- 
ing sees a band of intensity between T and M', as well 
as T and K' . This is absent in the singlet case, which 
can distinguish between these states. Such a band is also 
seen in the "singlet + triplet" J7(l)[7r, 0] state, which oth- 
erwise is again similar to the corresponding singlet case. 

The Z2 (S+T) spin liquids are shown in Fig. 12. 
Most of the structure factors closely resemble their singlet 
counterparts with important general differences. Here, 
we shall limit ourselves to pointing out these differences 
only. Generally, due to the triplet decoupling channels, 
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FIG. 10. (Color online). Dynamical structure factor for the (a) Z 2 [0,0]A, (b) Z 2 [0,0]£, (c) Z 2 [0,0]D, (d) Z 2 [n,n]B (e) 
Z2[0,7r]/3, (f) Z2[0,7r]a, (g) Z2[0,7r]7, and (h) Z2[7r,0]B singlet spin liquid ansatze. These structure factor are plotted from the 
centre of the extended Brillouin zone of the Kagome lattice, T, to the edge M\ to the corner K\ and back to the center T. 
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FIG. 11. (Color online). Difference between the dynamical structure factor for the (a) £7(1) [0,0], (b) C7(l)[0, 7r], (c) C/(l)[7r,7r], 
and (d) Z7(l)[7r, 0] "singlet + triplet" spin liquid ansatze. These structure factor are plotted from the centre of the extended 
Brillouin zone of the Kagome lattice, T, to the edge M\ to the corner K' , and back to the center Y. We note that since the 
"singlet + triplet" C/(1)[0, 7r] spin liquid state is gapped, it is unstable to instanton tunnelling events. 64 
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the spinon bands are spin-split. Hence, the sharp dispers- 
ing structures are somewhat lost and the scattering inten- 
sity becomes more diffuse compared to the singlet only 
ansatze, with increasing strength of the triplet terms. As 
with the £7(1) states, the intensity of low energy scat- 
tering does not go to zero at q = 0. Hence there is a 
slight enhancement of scattering at the EBZ centre com- 
pared to the singlet states. Also, for most of the states 
where there is a spin-gap in the structure factor, the gap 
magnitude becomes increasingly independent of the mo- 
mentum. 

Finally, we end this section by noting that the spin- 
liquid state Z2[0,7r]/3 seems to be the best candidate in 
the context of present inelastic neutron scattering exper- 
iments and DMRG calculations. We shall discuss this in 
some detail, along with the signatures of the other spin 
liquids, in the next two sub-section in the context of in- 
elastic neutron scattering and ESR experiments. 



E. Implications for Inelastic Neutron Scattering 
Experiments 

It is useful to compare the spin-structure factor ob- 
tained above with the results of the recent inelastic 
neutron scattering experiments on single crystals of 
Herbertsmithite. 9 These experiments show almost uni- 
form diffuse scattering over a large energy window. Fur- 
thermore, no spin-gap is observed down to 1.5 meV. It 
also shows no obvious signature of broken lattice symme- 
tries. 

With parameters from the mean-field solution, (as dis- 
cussed later in Appendix C) excitations are found for 
uj < 0.35 J. In contrast, the inelastic neutron scattering 
shows strong diffuse intensity up to the largest measured 
energy 11 meV ~ 0.65 J. However, the numerical values 
from our mean- field calculation can, at best, serve as a 
consistency check. 

Comparing with our calculated dynamical spin- 
structure factor with the above features of the exper- 
iment, we may infer that the £7(1) [71", 7r], £7(1) [71", 0], 
Z2[7r,7r]£?, and Z 2 [7r,0]£> states clearly appear to be in- 
consistent with the neutron scattering data. Further 
comparison of their ESR line shapes and peak distri- 
bution (see next section) may confirm/ invalidate this 
conclusion. 

In contrast, several of the [0,0] and [0, it] (both £7(1) 
and Z2) QSLs do have mostly- featureless diffuse scatter- 
ing within some energy window [co> m in, w max ]. However, 
in all but a few cases this window is narrow, and these 
states have well-defined features at lower or higher ener- 
gies and hence are inconsistent with experiments. 

The four states, namely Z2[0,0]D, Z2[0,7r]a, Z2[0,7r]7 
and Z2[0,7r]/3, have dispersive features with a very weak 
intensity, leading to an almost diffuse structure factor. 
The dispersive features are broadened and are in rel- 
atively poor contrast with the generally diffuse back- 
ground. However, we see some modulation in intensity 



for Z 2 [0,7r]7 state, and slightly better-defined disper- 
sive features in the Z2[0,0]D state. However, only the 
Z2[0,7r]/3 state shows evidence of a spin gap - Z2[0,7r]a 
is gapless. As already noted, the inclusion of the spin- 
rotation symmetry breaking perturbations makes spin- 
gap increasingly momentum independent in Z<$$,ii\f}. 
Further these perturbations also enhances the scatter- 
ing near the Brillouin zone centre somewhat in almost 
all cases (refer to Fig. 12). Noting that the DMRG 
calculations indeed stabilize a gapped QSL state (most 
likely Z2), in light of the inelastic neutron scattering mea- 
surements, it is tempting to suggest the Z2[0,7r]/3 state 
as a consistent choice for a candidate ground state for 
Herbertsmithite. However, we should note that the is- 
sue of existence of spin-gap is still not clear in Herbert- 
smithite experiments, and except for the lack of spin-gap, 
Z 2 [0,7r]a is also a possible candidate. 

Next, we describe the ESR absorption spectra for the 
QSLs that break spin-rotation symmetry. 



F. ESR Absorption Intensity for different 
spin-rotation symmetry breaking QSLs 

To observe the ESR spectra we couple the system with 
a Zeeman field via Hz z-Si. The ESR absorption in- 
tensity given in (25). We note that I(uj) oc uj at low 
temperatures, and peaks occur around the value of the 
Zeeman field, {uj = Hz) so absorption is most easily seen 
at the large values of uj and Hz- We align the Zeeman 
field along the z-axis, the axis along which the spin- 
rotation is broken. For energy scales uj near Hz-, the 
Zeeman term breaks SU{2) symmetry, with Uij along 
with the triplet terms. For this axis of the magnetic 
field, I xx {uj) = I yy (uj) due to the remaining £7(1) spin- 
rotational symmetry around the z-axis. 

We now discuss the ESR absorption intensity for the 
four U(l) and eight Z2 QSL in presence of the triplet 
channels. We focus on the small-intensity region to show 
contributions from the triplet terms; these occur as satel- 
lite peaks with smaller intensities, while the peak at 
uj = Hz has a much larger intensity. These spectra are 
shown in Fig. 13. 

Both the £7(1) [0,0] and £7(1) [0, n] spin liquids have an 
additional 'satellite' peak on either side of the main peak 
at uj = Hz- The peaks of the Z7(l)[0,7r] state have both 
higher intensity and are spaced more closely compared 
to the f7(l)[0,0] state. 

However, the Z7(l)[7r,7r] and Z7(l)[7r,0] states have 
wildly different ESR absorption spectra. There is absorp- 
tion over a broad range of Hz—uJ, with many subsequent 
satellite peaks almost forming a continuum. The inten- 
sity of the main peak is significantly diminished com- 
pared to the /7(1)[0,0] and £7(1) [0, 7r] spin liquids. As 
Hz — uj increases from zero, in the £7(l)[7r,7r] state, the 
primary satellite peaks have a monotonically decreasing 
intensity, while the £7(l)[7r,0] state features a single pair 
of prominent satellite peaks at a finite distance away from 
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FIG. 12. (Color online). Difference between the dynamical structure factor for the (a) Z 2 [0,0]A, (b) Z 2 [0,0]B, (c) Z 2 [0,0]D, 
(d) Z2[tv,7t]B (e) Z2[0,7r]y5, (f) Z2[0,7r]a, (g) Z2[0,7r]7, and (h) Z2[tv,0]B "singlet + triplet" spin liquid ansatze. These structure 
factor are plotted from the centre of the extended Brillouin zone of the Kagome lattice, T, to the edge M', to the corner K' , 
and back to the center T. 
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FIG. 13. (Color online). Electron spin resonance absorption for (a) 17(1) [0,0], (b) C7(l)[0,7r], (c) £7(1) [tt, tt], and (d) U(1)[tt,0] 
"singlet + triplet" spin liquid ansatze, as a function of the Zeeman field strength Hz and microwave field frequency oj. The 
viewpoint is along the Hz — oj axis. We note that since the "singlet + triplet" £7(l)[0,7r] spin liquid state is gapped, it is 
unstable to instanton tunnelling events. 64 
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zero. 

The Z2 spin liquids generally lead to a broadening 
of the absorption peaks. These spectra are shown in 
Fig. 14. The Z2[0,0]A state shows slight broadening 
of the all peaks around the base, with diminished inten- 
sity compared to the the U(l) case. The Z2[0,0]D state 
shows a similar response, though the satellite peaks are 
higher compared to the broadening of the main peak. 
The Z 2 [0,0]I? state has very broad absorption around 
the now significantly smaller satellite peaks. The largest 
satellite peaks of the Z2[tt,7t]B are broadened at low ab- 
sorption, but the intensity around the main peak drops 
dramatically and is almost zero. 

The Z2[0,7r]/3 state shows significant broadening of 
the main peak around the base, and nearly complete 
reduction of satellite peak intensity. The Z2[0,7r]a 
and Z 2 [0,7r]7 states have smaller amount of broadening 
around the main peak, where the satellite peaks are di- 
minished in intensity, but still visible. The intensity of 
the main peak of the Z2[ir,0]B is almost zero and the 
satellite peaks are replaced by small, broad domes of ab- 
sorption. 

The shape and structure of the satellite peaks provide 
another qualitative clue to distinguish between the spin 
liquid ansatze considered here, and immediate evidence 
of spin-rotation symmetry breaking in the kagome planes. 

Current ESR measurements down to as low as 5 K, 
however, find that the absorption intensity is dominated 
by the impurity spins below 20 K. 13 These Cu 2+ ions 
contribute in a nearly-paramagnetic fashion, and the in- 
tensity displays a Curie-like response ex 1/T. The line 
shape remains broad. In the regime where these impu- 
rity spins display a paramagnetic response, isolating the 
kagome-layer ESR response may prove difficult. An in- 
crease in sample purity can curtail this effect. However, 
at very low temperatures, these impurity spins may no 
longer behave in a paramagnetic fashion, and interact 
with the kagome planes with a strength ~ 10 K. 13,65 
One needs to consider the effect of the coupling of these 
impurities to the kagome spin liquid, and the resultant 
ESR spectrum. Such analysis is beyond the scope of this 
present paper, so we present here the intrinsic absorption 
of the spin liquid layers. 



VI. DISCUSSION 

We now summarize our results. In this work, we have 
attempted to address two important questions to account 
for the unusual phenomenology of the non-magnetic 
ground state of Herbertsmithite. Firstly, we have cal- 
culated the dynamical spin-structure factor for four U(l) 
and eight Z2 symmetric spin liquids derived from the 
former for the NN and NNN antiferromagnetic Heisen- 
berg model on an isotropic Kagome lattice. These spin 
liquids are allowed by spin-rotation, time-reversal and 
lattice symmetries of a Kagome lattice. We then con- 
sider the effect of small spin-rotation symmetry breaking 



perturbations (DM and Ising anisotropy) on the above 
spin liquids, and on the corresponding spin-structure 
factor. Furthermore, we calculate the ESR absorption 
spectra, which show non-trivial structures in presence of 
spin-rotation symmetry breaking. Since recent numeri- 
cal studies suggest that the ground state of strictly NN 
antiferromagnetic Heisenberg model is very sensitive to 
small perturbations of second neighbour exchange and 
Dzyaloshinsky-Moriya interactions, we expect that the 
the above perturbations may have important effects in 
the low energy features of the experimentally measurable 
spin-structure factor. 

Indeed, we find that the addition of the perturbations 
make the structure factor largely diffuse over an extended 
energy window throughout the Brillouin zone. Simi- 
lar scattering has been observed in the recent inelastic 
neutron scattering experiment on Herbertsmithite. We 
particularly find two Z2 spin liquid states, the so-called 
Z2[0,7r]/3 and Z2[0,7r]a states, whose spin-structure fac- 
tor features are qualitatively in conformity with the ex- 
periments. Only the Z 2 [0,7r]/3 state, however, exhibits a 
spin-gap in the structure factor. In the presence of the 
perturbations, this gap becomes increasingly momentum- 
independent. Noting that a similar gapped QSL was 
obtained as the ground-state of the NN antiferromag- 
netic Heisenberg Hamiltonian in recent DMRG calcula- 
tions, it is tempting to infer that the ^[0, tt]/3 state may 
be adiabatically connected to the ground state obtained 
in DMRG and also the ground state of Herbertsmithite. 
This prediction can be checked further by measuring the 
ESR absorption spectrum. Our present calculations sug- 
gest that the the Z 2 [0,7r]/3 phase shows a characteristic 
broadening of lines in ESR absorption spectra. 
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Appendix A: U(l) Spin Liquid Ansatze 

The flux patterns for the nearest-neighbour £7(1) spin 
liquids under consideration are shown in Figure 15. They 
consist of positive and negative real hopping terms that 
retain the translational symmetry of the Kagome lat- 
tice (17(1) [0,0] and 17(1) [tt, tt]) or break it (J7(1)[0,tt] and 
/7(1)M]). 
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FIG. 14. (Color online). Electron spin resonance absorption for (a) Z 2 [0,0]A, (b) Z 2 [0,0]B, (c) Z 2 [0,0]D, (d) Z 2 [n,n]B (e) 
Z2[0, 7r]/3, (f) Z 2 [0,n]a, (g) Z 2 [0,n]j, and (h) Z2[7r,0]B "singlet + triplet" spin liquid ansatze, as a function of the Zeeman field 
strength Hz and microwave field frequency oo. The viewpoint is along the Hz = oo axis. 
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FIG. 15. (Color online). Flux pattern for the real hopping 
terms of the (a) £7(1) [0,0], (b) C7(1)[0,tt], (c) £7(1) [tt, tt], and 
(d) C7(l)[7r, 0] singlet spin liquid ansatze. Hoppings are either 
positive (dark blue, solid) or negative (light brown, dashed) 
on nearest-neighbour bonds. 



Appendix B: Projective Symmetry Group 

Determining the allowed PSGs {GsS}, and thereafter 
the ansatze T~Lq, comes from constraints offered by the 
group multiplication table of the symmetry group SG. 
With group multiplication rules like AB = C that can be 
rewritten as ABC~ X = I, we can place constraints on the 
corresponding expression (GaA)(GbB)(GcC)~ 1 . First, 
there is no net physical transformation, so the expression 
must reduce to a gauge transformation. Second, each of 
the operations leaves 7~Lq invariant, and so does the final 
expression. Thus, (G A A)(G B B)(GcC)- 1 e IGG. 

As derived by Lu et. al, Gs are given in the singlet 
case by the following: 10 



G Tl {x,y,s) 


= v v 12 i 


(Bl) 


G T2 (x,y,s) 


= i 


(B2) 


G<r(x,y,s) 


= V12 9cr(s) 


(B3) 


Gc 6 (x,y,s) 


= vT x(x+1)/2 9c 6 (s'),s'G{u,v} 


(B4) 


G C6 (x,y,w) 


= vT X+y+XiX+1)/2 9C 6 (w) 


(B5) 


Gt 


= ir\ 


(B6) 



where we specify the lattice site i by the position of the 
unit cell to which is belongs R = xa. + yh, and the sub- 
lattice index s = u, v, w, as indicated in Fig. 2. 

The PSG can also be used to generate terms on 
symmetry-related bonds of an ansatze. The action of 
a lattice transformation S and an SU (2) gauge transfor- 
mation W on the quadratic terms of bond, given by Uij, 
is 



SUij -> C/s-i^s-iq) 



WUi. 



W(i)UijW\j). 



(B7) 



If GsS leaves the Uij invariant, as in the singlet case, 
then we can combine the above relations to generate Uij 
on symmetry-related bonds: 



G s (S(WijGUsU))- 



We will cover the case of the "singlet + triplet" PSG 
next. 



1. Projective Symmetry Group for "Singlet + 
Triplet" Ansatze 



Now, we discuss how time reversal and (improper) ro- 
tations affect "singlet + triplet" Uij, as in (15). We will 
see that the PSG {GsS} are the same as in the singlet 
case, and that the ansatze allow imaginary triplet terms 
as well. 

We will begin by determining how the standard 
fermionic time-reversal operator Tf changes the ansatz 
Uij, to properly capture the effect of time-reversal on the 
triplet terms. Tf = #K, where K is the complex conju- 
gation operator and 





(B9) 



Consider the action on a generic U^: 

I A R\ I - 

U = 




One can instead define a modified time-reversal operator, 
with an additional ir 2 gauge transformation acting on 



Uij, T : 



iT*L f , giving 




-U 



Since the C3 rotations are performed around the z-axis 
in both real space and spin space, neither the singlet nor 
the E Z ,Y Z triplet terms in our ansatze are affected by 
them. 

As all ansatz will projectively obey time-reversal sym- 
metry, we can use the form of Gt to simplify the ansatz 
before considering the effect of reflection. The action of 
T upon U^ is —U^ in both singlet and triplet cases, 
thus they affect the ansatze in the same manner. Thus, 
T 2 = +1 acting on the mean-field states, as before, and 
we have the same condition ±Gt(^) 2 £ IGG. 

The action in spin space of the reflection, cr, is to 
flip the sign of our triplet components, which are the 
z-components of a pseudo-vector. Since a gauge trans- 
formation, on physical grounds, can only mix singlet and 
triplet terms among themselves, this action commutes 
with the SU (2) gauge transformations. Only the lattice 
part of cr enters in the commutation relation between cr 
and the gauge transformations Gs- Solving for the time- 
reversal gauge transformation Gt as in the singlet case, 
the non-zero mean-field anatze will all have the same 
Gt = ir 1 . 10 For an ansatz to be time-reversal invariant 
we must have 



(B8) 



-U t . 
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With this restriction, Uij can be parametrized by Uij = 
72T 2 + 73T 3 , where 72,3 £ C. We note that the real com- 
ponents of 72,3 are coefficients of the singlet pairing and 
hopping, respectively, while the imaginary components 
are coefficients of triplet pairing and hopping. Reversal 
of the sign of the triplet coefficients can be performed by 
taking Uij -> U\- = Uji. 

With this, we have the same projective symmetry 
group that was derived for singlet spin liquids with the 
symmetry of the Kagome lattice. 10 However, the ansatz 
that can be generated will differ, due to the effect of a 
upon the triplet terms. 

We may generate all nearest and next-nearest- 
neighbour bonds from just one by using the symmetry 
operations of the Kagome lattice. For each, we may im- 
pose a constraint on the allowed ansatz for these bonds 
by devising a non-trivial lattice symmetry operation that 
takes Uij back to itself. With the effect of lattice trans- 
formations (B8) along with the transformation under cr, 
we can derive the following constraint on the nearest- 
neighbour Uij\ 

gai^gcei^gcsMU^glj^gl^gljiv) = C/ NNj 

(BIO) 

and a similar constraint on the next-nearest-neighbour 



ga(u)gc 6 ( u )UNNNgc 6 ( v )9a( w ) = ^NNN- (Bll) 

The particular bonds for Z7nn and Z7nnn are shown in 
Fig. 1. For a given Z2 ansatz in Table III, these may 
restrict the structure of Uij to disallow hopping (r 3 ) or 
pairing (r 2 ) terms, but does not place any conditions on 
whether these terms are real or imaginary. This differs 
from the singlet case, where the effect of a upon the 
triplet terms was not considered, and the ansatze become 
restricted to only singlet terms. Thus, spin-rotational 
symmetry-breaking triplet terms are allowed within the 
same PSG for the Kagome lattice. Since U^ is now no 
longer Hermitian, Uji 7^ Uij, and the direction of the 
bonds matter. They are shown in Fig. 1. 



SL Label 


7712 


gc 6 


go- 


Z 2 [0, 0]A 


+1 


u,v,w : I 


U, V, w 


: I 




-1 


u,v,w : I 


U, V, w 


: I 


Z 2 [0,0]B 


+1 


u,v,w : ir 3 


U, V, w 


: I 


Z 2 [0,7r}a 


-1 


UjVjW : ir 3 


U, V, w 


: / 


Z 2 [0,0]D 


+1 


u,v,w : ir 3 


u,v,w : 


IT 3 


Z 2 [0,7r] 7 


-1 


u,v,w : ir 3 


u,v,w : 


IT 3 


Z 2 [tt,tt]B 


+1 


u,v : I w : ir 1 


u,v,w : 


IT 3 


Z 2 [k^]B 


-1 


u,v : I w : ir 1 


u,v,w : 


IT 3 



TABLE III. Parameters 7712, #c 6 and gcr characterizing the 
PSG {GsS} of spin liquid (SL) states, as given in (B6). 



Appendix C: Mean-Field Results 

As mentioned in Sec. II, we consider the mean-field so- 
lution for each of the Z2 spin liquid ansatze, with nearest- 
neighbour perturbations D, A and the next-nearest- 
neighbour J 2 . While we do not expect the self-consistent 
mean-field theory to give a quantitatively correct phase 
diagram, we can still gain some insight from the results. 
In particular, we would like to understand the represen- 
tative mean- field parameter values U^ when the states 
in Table II are stabilized. 

We focus on the parameter regime D / J, A/ J £ [0,0.5). 
Furthermore, we will fix cu_ as the overall nearest- 
neighbour energy scale, taking the regime J2/00- £ [0, 2). 
In this way, we separate the effects of D and A from J2 
as much as possible. We begin by considering general 
trends across all ansatze. 



1. General Results 

DM interactions: Here, we consider the contributions 
from the x+ and 77+ channels in Eq. (18), in comparison 
to contributions from X- and 77 _ channels. The self- 
consistent theory suggests that for U(l) and Z2 states 
labelled as [0, 0] and [0, n] states, the \- and 77- channels 
have the dominant contributions, and the ratio of triplet 
to singlet terms is a nearly linear function of D / J within 
our parameter window < D, A < J/2. The [7r,7r] and 
[-7T, 0] states, however, have stable mean-field states where 
X+ and 77+ are significant, and dramatically small triplet 
values, particularly for D/J < 0.2. 

Ising Interaction A: This term has a negligible effect 
on the Z2[0,0] and Z2[0,7r] states, since x+ and 77+ are 
not relevant. For Z2[7r,7r] and Z2[7r,0] spin liquids, how- 
ever, small values of A/ J make x+ and 77+ channels less 
relevant. Since uu+ = (J — A)/8 + O(D), increasing A 
from zero will actually decrease so these channels 
make little contribution to Hq in the mean-field theory. 
Within the mean-field theory, these small values of A 
yield no qualitatively new behaviour. 

Z2 Spin Liquids and Next- Nearest- Neighbour Terms: 
For the six (all except Z2[0,0]£? and Z2[0,7r]a) Z2 spin 
liquids in Table II that are stabilized only in presence 
of next- nearest- neighbor terms, within mean-field the- 
ory, we must have 0.5 < J2/0J- < 1.2. The ^[0,0] and 
Z2[0,7r] states have a next-nearest-neighbour \2 hopping 
parameter that does not stabilize the Z 2 spin liquids from 
their U(l) parents. Generally, |x2 1 increases with J2 in a 
monotonic and nearly-linear fashion. However, at a large 
value of J2 ~ 1.2co>_, the Z2[0,7r]7 state is stabilized with 
small |xi|, small 1 772 1 , and large |x2 1 - We expect that such 
a state is unlikely to be realized in models where J2/J 
is small. The Z2[tt,7t]B state is stabilized with a similar 
jump at J 2 ~ 0.9o;_ with a large \\2\ value, and a small 

Out of the other four, in three Z2 spin liquids 
(Z2[0,0]A, Z2[0,0]D and ^[0, 7r]/3), 77 and \2 increases 
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FIG. 16. Heuristic mean- field phase boundaries in the J2- 
D plane of ^[0, 0] states, where J2 is next-nearest-neighbour 
Heisenberg coupling, D the nearest-neighbour DM interaction 
in the z direction, and cj_ is defined in (18). 
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FIG. 17. Heuristic mean- field phase boundaries in the J2- 
D plane of Z2[0,7r] states, where J2 is next-nearest-neighbour 
Heisenberg coupling, D the nearest-neighbour DM interaction 
in the z direction, and uj- is defined in (18). 



within mean- field theory. The Z 2 [{),{)\B phase is stabi- 
lized throughout the phase diagram, while the Z 2 [0,7r]a 
state is never stabilized within mean field theory of the 
current Hamiltonian. The other Z 2 spin liquids are sta- 
bilized for sufficiently large J 2 . 

Fig. 16 shows the U (1)-Z 2 phase boundary in the J 2 -D 
plane for the Z 2 [0, 0] spin liquids, Fig. 17 for the Z 2 [0, n] 
spin liquids, and Fig. 18 for the Z 2 [7r,7r]B and Z 2 [k,0]B 
spin liquids. The next-nearest-neighbour interaction J 2 
stabilizes all other six Z 2 spin liquids. We see that the 
Z 2 [0,0U, Z 2 [0,7rl/3 and Z 2 \tt,0]B phases are stabilized 
0.5 1 1 1 1 . 1 1 

\ Z 2 [n,0]B 

\Z 2 [n,n]B 




FIG. 18. Heuristic mean- field phase boundaries in the 
J2-D plane of Z^^\B and Z2[tt,0]B states, where J2 is 
next-nearest-neighbour Heisenberg coupling, D the nearest- 
neighbour DM interaction in the z direction, and uj- is defined 
in (18). 



for smaller J 2 than the others, while being destabilized by 
the Dzyaloshinsky-Moriya interaction D. However, the 
Z 2 [0,0]D, Z2[0,7r]7 and Z 2 [tt^7t]B phases are stabilized 
by D, despite requiring larger J 2 . 



monotonically with J 2 beyond a critical value of J 2 /J 
which itself depends on the type of spin liquid in con- 
sideration. While the Z 2 [tt,0]B sees this behaviour for 
smaller values of J 2l it also undergoes a jump around 
J 2 ~ 0.9cj_, similar to the Z 2 \k,t\\B state. 



2. Mean-Field Phase Boundaries 

Here we present the phase boundaries in the J 2 -D 
plane where these Z 2 spin liquid phases can be stabilized 



Appendix D: Dynamical Structure Factor of 
Singlet +Triplet Ansatze 

In the main text, we showed the changes in the dy- 
namical structure factor from the singlet ansatze to 
when small triplet terms are added. Here we present 
the trace of the dynamical structure factor matrix, i.e., 
(S xx + S yy + S zz )/3, for these triplet ansatze. We show 
the dynamical structure factor for the U(l) and Z 2 spin 
liquids in Fig. 19. 
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